Mon. Not. R. Astron. Soc. 000,[l]|2l(2005) Printed 5 February 2008 (MN MbX style file v2.2) 



X-rays from the radio-quiet quasar PG 1407+265: relativistic jet or 
accretion disc emission? 

L. C. Gallo 

Max-Planck-Institut fur extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany 



Accepted. Received. 



ABSTRACT 

We present two XMM-Newton observations of the luminous (L^ > lO^^ergs"^ ), radio- 
quiet quasar, PG 1407+265, separated by eleven months. The data indicate two distinct states: 
a highly variable, bright state (first epoch); and a quiescent, low-flux one (second epoch). 
During the low-flux state the spectrum is consistent with a single, unabsorbed power law. 
However, during the brighter state a highly variable, steep component is statistically required. 
Contemporaneous UV data from the Optical Monitor allow an estimate of the optical-to-X-ray 
spectral index (aox), which appears typical of radio-quiet quasars during the low-flux state, 
but extremely flat during the high-flux state. The XMM-Newton data can be described as origi- 
nating from a combination of jet and accretion disc processes, in which the (relativistic) X-ray 
jet only works intermittently. The scenario could help describe some of the complexities seen 
in the broadband spectral energy distribution of PG 1407+265, such as weak high-ionisation 
emission lines, strong Fe II, unbeamed continuum, and the weak radio emission relative to the 
optical. 

Key words: galaxies: active - galaxies: jets - galaxies: nuclei - quasars: individual: 
PG 1407+265 - X-ray: galaxies 



1 INTRODUCTION 

The radio-quiet quasar, PG 1407+265 (z — 0.94), is best known 
for constant display of enigmatic behaviour. A detailed analysis of 
the spectral energy distribution (SED; McDowell et al. 1995) re- 
vealed that PG 1407+265 exhibited a normal non-variable, radio- 
quiet continuum, but weak and blueshifted high-ionisation emis- 
sion lines. In contrast, the UV and optical Fe II appeared rather 
strong. A hand-full of weak-line quasars have since been discov- 
ered (e.g. Leighly et al. 2004; Hall et al. 2004; Reimers et al. 2005), 
but no general consensus has been reached as to the nature of the 
physical mechanisms involved. 

A relativistically beamed continuum, as in a blazar, could ex- 
plain the absence of emission lines. For PG 1407+265 this expla- 
nation is normally dismissed as the SED of the quasar is strikingly 
similar to that of radio-quiet AGN (Elvis et al. 1994). Moreover, 
PG 1407+265 is relatively radio weak with an optical-to-radio flux 
ratio of 3.43 (Kellermann et al. 1989) and, by definition, consis- 
tent with radio-quiet AGN. However, a multi-year radio observing 
campaign of PG 1407+265 with the VLBA revealed strong evi- 
dence of a radio jet with a highly relativistic speed (Dopper factor 
> 10) and likely viewed within a few degrees of pole-on (Blundell 
et al. 2003). Blundell et al. suggested that the pole-on orientation 
could explain some of the opticalAJV properties, assuming that the 
line-of-sight was into the broad line region, thus the emission line 
spectrum is diluted by the accretion disc continuum. 



Although PG 1407+265 is individually curious, a better under- 
standing of it could have broader implications. The debate whether 
the radio characteristics of AGN show a bimodal or continuous 
distribution is ongoing (e.g. Goldschmidt et al. 1999; White et al. 
2000; Ivezic et al. 2002; Cirasuolo et al. 2003). The discovery of 
objects which possess characteristics of both groups are obviously 
of interest. It is clear that all AGN are capable of producing ra- 
dio emission, hence able to generate jets on some level. Observa- 
tional support for this comes from sensitive, high-resolution radio 
imaging of radio-quiet AGN, evidencing jets in some objects (e.g. 
Blundell & Beasley 1998; Ulvestad et al. 2005). The simplest con- 
clusion is that the radio mechanism in radio-quiet and radio-loud 
AGN are identical, except that the bulk kinetic power is smaller in 
radio-quiet objects. 

To some extent, these observations form the bases of the 
'aborted jet' model which has become fashionable in AGN circles. 
The idea being that short or aborted jets are produced in radio- 
quiet AGN and the high-energy jet particles, which are expelled 
close to the black hole environment, illuminate the accretion disc 
(e.g. Henri & Petrucci 1997; Ghisellini et al. 2004). The presence 
of these energetic particles can account, in whole or in part, for the 
high-energy continuum traditionally associated with the accretion 
disc corona. 

PG 1407+265 could evolve to be an important laboratory in 
investigation of the radio-quiet/loud dichotomy (or lack of it), be- 
cause the quasar clearly possesses characteristics of both classes: 
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it is an example of a radio-quiet AGN which has successfully 
launched a relativistic radio jet to large distances. In this paper, we 
examine the X-ray properties of this quasar to study the behaviour 
of possible jet emission closer to the central engine. 



2 OBSERVATIONS, DATA REDUCTION, AND FITTING 
PREPARATIONS 

PG 1407-1-265 was observed wi\hXMM-Newton (Jansen et al. 2001) 
on 2001 January 23 (revolution 206) and 2001 December 22 (rev- 
olution 373) for approximately 71 and 42 ks, respectively. At both 
epochs, all on-board instruments functioned normally. The EPIC 
pn (Striider et al. 2001) and MOS (MOSl and M0S2; Turner et 
al. 2001) cameras were operated in full-frame mode with the thin 
filter in place. The Reflection Grating Spectrometers (RGSl and 
RGS2; den Herder et al. 2001) also collected data during this time, 
as did the Optical Monitor (OM; Mason et al. 2001). The RGS data 
were presented in an examination of line-of-sight absorption to- 
ward PG 1407-1-265 (Fang et al. 2005), and will not be discussed 
here. 

The Observation Data Files (ODFs) were processed to pro- 
duce calibrated event lists using the XMM-Newton Science Analy- 
sis System (XMM-SAS v6 . 4 . 0) and the most recent calibration 
files. Due to anomalies in the Mission-Elapsed-Time during the 
first observation, the default values of the the environment variable 
SAS_OBTjyiET_FIT^ were changed accordingly to "0,1,0.5,0.9,1". 
Otherwise, standard processing practices were followed. Unwanted 
hot, dead, or flickering pixels were removed as were events due to 
electronic noise. Event energies were corrected for charge- transfer 
losses, and time-dependent EPIC response matrices were generated 
using the SAS tasks ARFGEN and RMFGEN. Light curves were ex- 
tracted from these event lists to search for periods of high back- 
ground flaring. Some flaring was evident during revolution 206, 
consequently the data during these intervals were simply ignored. 
The source plus background photons were extracted from a circular 
region with a radius of 35", and the background was selected from 
an off-source region with a radius of 50" and appropriately scaled 
to the source region. Single and double events were selected for the 
pn spectra, and single-quadruple events were selected for the MOS. 
Pile-up was examined at both epochs and deemed negligible in all 
instruments. 

The OM worked in imaging mode at both epochs. In total, 
eight images were taken during revolution 206, and ten during rev- 
olution 373. Only the UVW2 filter (1800 - 2250 A) was used 
with typical exposures between 3280 — 5000 s. A log of the XMM- 
Newton observations is provided in TableQ 

In the following analysis, the spectra are grouped such that 
each bin contains at least 20 counts, and fitting is performed 
using XSPEC V12.2.0 (Amaud 1996). Fit parameters are re- 
ported in the rest frame of the object and quoted errors corre- 
spond to a 90% confidence level for one interesting parameter 
(i.e. Ax^ = 2.7 criterion). K-corrected luminosities are derived 
assuming isotropic emission, a value for the Hubble constant of 
Ho=70 km Mpc~^, and a standard cosmology with Q,m = 0.3 
and = 0.7. 

A value for the Galactic column density toward PG 1407-1-265 
of 1.38x10^'' cm~2 (Elvis etal. 1989) is adopted in all of the spec- 
tral fits. Initial spectral fits were performed with a free absorption 

^ http://xmm.vilspa.esa.es/sas/current/doc/idxEnvironments.html 



Table 1. Log of XMM-Newton observations of PG 1407+265. The XMM- 
Newton revolution number of when the observation was conducted is given 
in column (1) and the observation id in column (2). The instrument is shown 
in column (3). The total amount of useful exposure (GTI) is reported in col- 
umn (4), and the estimated number of source counts is reported in column 
(5). The combined exposure for all images is given for the OM. The energy 
bands considered are 0.25 — 10 kcV and 0.3—10 kcV for the pn and MOS, 
respectively. 



(1) 


(2) 


(3) 


(4) 


(5) 


Rev. 


Obs. ID 


Instrument 


Exposure ( s) 


Counts 


206 


0092850101 


pn 


51458 


125750 






MOSl 


62925 


35331 






MOS2 


62821 


35419 






OM 


36000 




373 


0092850501 


pn 


35098 


30980 






MOSl 


41017 


8979 






MOS2 


41007 


8784 






OM 


38800 





parameter set to the redshift of the source. In determining that the 
intrinsic absorption was imperceptible (< 5 x 10^^ cm"'^ ), subse- 
quent spectral modelling included only the Galactic component. 

Recent work on the pn calibration show that the spectral re- 
sponse is reasonable down to 0.25 keV (F. Haberl priv. comm.). We 
confirmed that inclusion of the 0.25 — 0.3 keV data did not hamper 
the fit; therefore we elected to utilise the lower energy range in our 
analysis. 

On the low-energy end, the MOS are limited to 0.3 keV to 
minimise cross-calibration uncertainties between the three cameras 
(Kirsch 2005). At high energies, all three spectra are restricted at 
10 keV where background contamination becomes significant. 

3 THE BRIGHT-STATE JANUARY 2001 OBSERVATION 
3.1 Spectral analysis 

Fitting the pn and MOS spectra separately revealed relatively good 
agreement within known cross-calibration uncertainties (Kirsch 
2005); therefore for the initial phenomenological fits all of the 
EPIC data were fitted together, while residuals from each instru- 
ment were examined separately to distinguish gross discrepancies. 
A moderately steep power law (F = 2.32 ± 0.02) fitted the ob- 
served 1.3 — lOkeV band (2.5 — 19.4 keV rest frame) relatively 
well (xi5/dof = 0.96/699). A search for a narrow (ct < 150 eV) 
iron emission line between 6 — 7keV yielded no detection, with 
the maximum equivalent width of any such feature being EW < 
12 eV. 

Extrapolation of the power law fit to 0.25 keV in the pn 
and 0.3 keV in the MOS degraded the fit quality (xS/dof = 
2.84/1047), and showed a gradual excess in the residuals toward 
lower energies (Figure 0. Several models were attempted to de- 
scribe the broadband (~ 0.49 — 19.4 keV rest frame) spectra of 
PG 1407-1-265 in a phenomenological manner (Table|2j. 

First, the single power law was formally fitted to the broad- 
band spectra. This was an improvement over the extrapolated, high- 
energy power law, but, predictably, the data at higher energies 
(where statistics are lower) were underestimated. 

Broken power law, double power law, and blackbody plus 
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Figure 1. The residuals remaining in tlie pn (top panel) and MOS (lower 
panel; MOSl red, filled circles; MOS2 green, open squares) data from 
revolution 206, when extrapolating the 1.3 — lOkoV band power law 
(r fS! 2.32) to lower energies. The data are binned for display purposes. 



Table 2. Combined fit to all of the EPIC data of PG 1407+265 during rev- 
olution 206. The model, quality of fit, and rest frame model parameters are 
given in columns (1), (2), and (3), respectively. 



(1) 

Model 


(2) 


(3) 
Parameters 


power law 


1.11/1047 


r = 2.40 ±0.01 


broken 
power law 


1.07/1045 


Ti = 2.41 ±0.01 
Eb = 4.43lg-^2 keV 
F, - 2 27+''-°'* 

1 2 — og 


blackbody 

plus 
power law 


1.09/1045 


kT = 204^11 eV 
r = 2.36 ±0.01 


double 
power law 


1.08/1045 


Ti = 2.47 ±0.02 
r2 = 1.75 ±0.18 



power law models, all provided comparable improvement over the 
single power law fit, and the inclusion of a second continuum com- 
ponent (two additional free parameters) was a significant improve- 
ment. For example, when tested with an F-test, the broken power 
law model resulted in an F-test probability of 10~^. We note that 
the double power law fit presented in Table|2|was not unique. Dif- 
ferent initial parameters resulted in slightly different final parame- 
ters, but comparable xt- 

Remarkable was the high energy where the apparent change 
in spectral slope occurs in all of the models. This is seen in the 
high break energy (Eb ~ 4.43 keV rest frame) in the broken power 
law, and also manifested in the high temperature of the blackbody 
component (kT ~ 204 eV). As in the broken power law fit, the 
steeper component dominates the spectrum up to rather high en- 
ergies in the double power law model, as well. It is unlikely that 
the high temperature or high break energy are physically meaning- 
ful, rather they instead indicate the presence of another component. 
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Figure 2. The 0.2- lOkeV pn fight curves of PG 1407+265 during the Jan- 
uary 2001 (rev. 206) observation (top curve) and December 2001 (rev. 373) 
observation (bottom curve). There is clear variability during the brighter 
observations, whereas the light curve is consistent with a constant during 
revolution 373. The data are binned in 500 s, and s on the time axis marks 
the start of each observation. 



which cannot be uniquely determine with the current statistics (see 
Sect.lO. 

Considering the broken power law fit, the observed 0.3 — 
lOkeV flux, corrected for Galactic extinction, was about 5 x 
10"^'^ ergcm"^ . The corresponding rest frame luminosities 
in the 0.3 - 2 and 2 - lOkeV range were 2.14 x 10'*'^, and 
8.85 X 10''"' ergs"' , respectively. The 2 — lOkeV luminosity was 
comparable to the 1993 measurement made from an ASCA obser- 
vation (L2-10 ~ 7.5 X 10**^, Reeves & Turner 2000). 



3.2 Variability 

3.2.1 Broadband X-ray flux variability 

Examining the observed 0.2 — 10 keV light curve in 500 s bins re- 
vealed significant variations on the order of ±20% about the mean 
count rate over the duration of the revolution 206 observation (see 
Figure |2}. The variability is even more impressive when we re- 
call that the average luminosity of PG 1407+265 was in excess of 
lO'*^ ergs"' , and that the quasar is radio-quiet. 

The most rapid rise occurred about 17ks into the observa- 
tions when the count rate increased about 10% in ~ 9 ks. Adopt- 
ing the blackbody plus power law fit resulted in the most con- 
servative estimate of the luminosity rate of change. Accordingly, 
the change in count rate corresponds to a luminosity change of 
AL « 2.84 X 10'*^ ergs"' in a quasar rest frame time interval of 
At « 4600 s. Following Fabian (1979; see also Brandt et al. 1999) 
the radiative efficiency of PG 1407+265 was rj ^ 0.29. Assuming 
a uniform, spherical emission region the measured rj is compara- 
ble to the efficiency expected from a maximally rotating Kerr black 
hole (77 ~ 0.3, Thorne 1974). The demand on black hole spin can 
be relaxed if asymmetric emission processes, such as relativistic 
Doppler boosting or light bending, play an important role in the 
energy output of PG 1407+265. 



4 L. C. Gallo 




0.2-0.7 keV counts s" 




0.2-10 keV counts s' 




-20000 -15000 -10000 



10000 15000 20000 



Delay (s) 



Figure 3. The count rate correlation between the 0.7 — 10 keV and 0.2 — 
0.7keV bands is shown in the top panel. The cross-correlation function 
between the two bands is shown in the bottom panel. The positive offset 
indicates that the hard band follows the soft band by about 3000 s. The data 
are binned in 500 s. 



3.2.2 Interband variability 

Initially, variability was examined in three sub-bands: 0.2 — 0.7, 
0.7 — 2, and 2—10 keV. It was established that there was no dis- 
cernible spectral variability or significant differences in the ampli- 
tude of the fluctuations between the 0.7 — 2keV and 2 — lOkeV 
bands, thus to improve statistics, the two bands were combined. 
Therefore, light curves are examined in the 0.2 — 0.7 keV (0.4 — 
1.35 keV rest frame) and 0.7- 10 keV (1.35-19.4 keV rest frame) 
bands. 

It was immediately noticeable that the degree of variability 
was much greater in the lower energy band. The fractional variabil- 
ity (JVar) in the 0.2 - 0.7 keV band was 7.33 ± 0.64%, whereas 
between 0.7 - lOkeV, Fvar = 3.34 ± 0.73% (uncertainties esti- 
mated following Edelson et al. 2002). 

The count rates measured in the two bands were well- 
correlated (Figure|3}. As such, we examined possible lags between 
the two bands by calculating the nominal cross-correlation func- 
tion. With the light curves in 500 s bins, there does appear to be a 



Figure 4. The hardness ratio (HR) as a function of 0.2 — lOkeV count 
rate during revolution 206. The spectrum becomes softer as the intensity 
increases. The pn data are binned in 500 s intervals. 



lag (Figure|3}, such that the low-energy band leads the hard energy 
band by ~ 3000 s. 

The significance of the lag should be treated with caution. 
When we cross-correlated the 0.2 — 0.7 keV band with 0.7 — 2 keV 
and 2 — 10 keV light curves separately, no lag was detected. It could 
simply be that the combined 0.7— 10 keV light curve produces suf- 
ficiently high signal-to-noise to make a detection, but nevertheless 
we do not consider the possible lag a robust detection. 



3.2.3 X-ray spectral variability 

As discussed in Section 13.2.21 there was clear spectral variability 
with the amplitude of the fluctuations being much larger in the 
0.2 - 0.7 keV band then in the 0.7 - lOkeV band. We investi- 
gated the spectral variability with time by calculating the hardness 
ratio (HR = H — S/H + S, where H and S are count rates in 
the hard and soft bands, respectively) variability curve. The hard- 
ness ratio curve was inconsistent with a constant (xl = 1.58) and 
correlated with the broadband count rate, such that the spectrum be- 
comes harder with diminishing intensity (Figure |4j. In general, as 
the flux gradually decreased over the ~ 70 ks observation, the spec- 
trum became harder. To investigate this short-term, flux-dependent 
spectral variability during revolution 206, we divided the observa- 
tion into half: the high-flux state (first half) and a low-flux state 
(second half). 

Initially, we fitted both pn spectra with a power law above 
3keV only. We found that this produced a good fit with a com- 
mon flux and photon index (F ~ 2.25), indicating that the spec- 
tral variability was dominated by the lower-energy component. Ex- 
trapolating this power law to lower energies, we determine that the 
high-flux spectrum appeared significantly steeper (Figure|5}- 

To investigate this further, we fitted both spectra simultane- 
ously with a broken power law. We maintained the break energy 
(Et) and high-energy power law (r2) common for both models 
and allowed only the low-energy power law to vary separately. This 
produced a good fit to the data (xt ~ 1-02) and illustrated the na- 
ture of the spectral variability (see fit parameters in Table |3}: the 
high-energy spectrum was consistent in flux and spectral shape; 
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Figure 5. The ratio from extrapolating a 3 — 10 keV power-law (F « 2.25) 
to lower energies for the high-flux (black crosses) and low-flux (red dots) 
pn spectra during revolution 206. The data are binned for display purposes. 



Table 3. A broken power-law fit to the high- and low-flux spectra of 
PG 1407-H265 during revolution 206. The energy band utilised is 0.25 — 
lOkeV (0.49 — 19.4 keV in the rest frame). The model parameter is 
shown in column (1), n is the model normalisation at IkeV in units of 
10~* photons keV~^ cm~^ s^^. The values during the high-flux and 
low-flux state are given in columns (2) and (3), respectively. Ej, and T2 
are fixed between the high and low state. 



(1) 


(2) 


(3) 


Model 


High-flux 


Low-flux 


Parameter 






Ti 


2.49 ±0.01 


2.39 ± 0.01 


Eb (keV) 


6.14+1-^4 


6.14+}-^^ 


r2 


9 97+0.09 
-0.11 


9 97+0.09 
^•^'-0.11 


n 


10.4 ±0.1 


9.7 ±0.1 



however, at low energies the variability seemed to be characterised 
by a change in shape as well as flux. 



3.2.4 Possible ultraviolet variability 

Eight images of PG 1407+265 were obtained in the UVW2 filter 
(rest frame 930 — 1160 A range) with the OM during the first 35 ks 
of the X-ray observation. For each image, a flux density was esti- 
mated from the count rate using the conversion factor derived by 
Chen (2004). The average UVW2 flux density during revolution 
206 was approximately 7.54 x 10~^^ ergs~^ cm~^ A 

We emphasise that within the flux uncertainties {^^ 5%) there 
were no significant variations in the UV light curve. Formally, the 
OM light curve was consistent with a constant. However, the appar- 
ent trend in the light curve was intriguing. Taking liberty to spec- 
ulate, when the UV light curve was overplotted (in time) on the 
X-ray light curve there appeared to be a common trend in the fluc- 
tuation (Figure|6}. The relationship between the curves could even 
be improved if we allowed for a slight time shift in one of them 
(e.g. if the UV lead the X-rays by a few kiloseconds). Admittedly, 



Figure 6. The UVW2 normalised light curve overplotted (in time) on a 
segment of the 0.2 — lOkcV normalised light curve during revolution 206. 
The left axis correspond to the UV data and the right axis to the X-rays. 
The s on the time axis marks the start of the X-ray observation. The pos- 
sible common trend could be strengthened by a slight time shift in the light 
curves. N.B. The figure should be considered illustrative, as the fluctuations 
in the UV are really much less then in the X-rays (note the values on the 
respective axes). 



Table 4. Combined fit to all of the EPIC data of PG 1407+265 during revo- 
lution 373. The energy band utiUsed is 0.25 - 10 keV (0.49 - 19.4 keV in 
the rest frame). The model, quality of fit, and rest frame model parameters 
are given in columns (1), (2), and (3), respectively. 



(1) 

Model 


(2) 
xS/dof 


(3) 
Parameters 


power law 


0.99/830 


r = 2.21 ±0.01 


broken 
power law 


0.98/828 


Fi - 2 36+"-^° 

1 1 — Z.OD_Q Q4 

= 1.30±0.19keV 

r2 = 2.19 ±0.02 


blackbody 

plus 
power law 


0.97/828 


kT = 108 ± 13 eV 
r = 2.19 ±0.02 


double 
power law 


0.98/828 


1 1 — Z.D(_g gg 

r, = 2.iil°:f. 



Figure|6|can appear deceiving, and we reiterate that the UV fluctu- 
ations are on the level of < 5%, whereas the X-rays variations are 
about 30%. Therefore, Figure |6| should be considered illustrative, 
in absence of higher quality UV data. 



4 THE LOW-FLUX DECEMBER 2001 OBSERVATION 
4.1 Spectral analysis 

A power law (F = 2.19 ± 0.04) provided an acceptable fit 
(xiJ/dof = 0.90/454 to the 1.3 - lOkeV (2.5 - 19.4 keV rest 
frame) EPIC spectra of PG 1407+265 during revolution 373. As 
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with the eariier observation, no narrow iron emission line was de- 
tected. 

Extrapolating the high-energy power law to lower energies 
was a good approximation to the broadband spectra. Indeed, fitting 
the single power law over this entire band did not result in a sig- 
nificant change of the photon index (AF ~ 0.02). The addition of 
a second continuum component did improve the fit over the single 
power law (F-test probability of ~ 10^"*), but not as significantly 
as in the revolution 206 data. 

Adopting the broken power law fit, resulted in an unabsorbed 
0.3 - lOkeV flux of about 2 x 10"^^ ergcm"^ s"^ . The corre- 
sponding rest frame luminosities in the 0.3 — 2 and 2 ~ lOkeV 
range were 7.29 x 10''^, and 3.95 x 10''^ ergs~^ , respectively. 

4.2 Timing analysis 

The timing behaviour of PG 1407-1-265 during revolution 373 
was much more in-line with the expected behaviour of a high- 
luminosity quasar. The 0.2 — lOkeV light curve (Figure |2j was 
completely consistent with a constant over the ~ 40 ks observa- 
tions (xS/dof = 0.81/76). Similarly, the UV light curve was also 
constant. 

The X-ray spectral variability was also negligible. The frac- 
tional variability in the two X-ray sub-bands (i.e. 0.2 — 0.7 keV 
and 0.7 — 10 keV) was on the 3% level, and the hardness ratio was 
also rather constant over time (x^ /dof — 1.11/76). 

5 COMPARISON BETWEEN THE TWO OBSERVATIONS 
5.1 Optical-to-X-ray spectral index 

The UVW2 flux density of PG 1407-1-265 diminished from (7.54± 
0.03) x 10"^^ erg s"^ cm"^ A"^ during revolution 206 to (6.82± 
0.05) x 10~^^ erg cm~^ A ^ during revolution 373. The frac- 
tional decrease in the AGN UV flux density was about 10% over 
the eleven month period. Over the same time, the average broad- 
band X-ray flux diminished by ~ 60%. 

A serendipitous, optically bright object (UVW2 ~ 13.51) 
was identified in the OM field-of-view in two images at each epoch. 
The object is possibly associated with the F8 star HD 124732, but 
this is not conclusive given the several arc minute discrepancy in 
positions. Of importance was that the flux of the object (2.01 x 
10"^"* erg cm~^ A ) was equal at both epochs (revolution 
206 and 373) to within about 1%. This clearly indicates that the 
different UV fluxes measured in PG 1407-1-265 arise from changes 
in the AGN and not from instrumental effects of the OM. 

The optical-to-X-ray spectral index, Oox, is the slope of 
a hypothetical power law extending between the optical (UV) 
and X-ray continuum in AGN: aox = log{fx/ fo) /log{ux/i'o)- 
In the standard definition fx and /o are the intrinsic flux den- 
sities at 2 keV and 2500 A, respectively, consequently a'^^ = 

0.384io^,(/...W/,5ool). 

Typical values of qJ,^ for unabsorbed, radio-quiet AGN are 
between —1.2 and —1.8 with the slope steepening with increasing 
UV luminosity (e.g. Strateva et al. 2005). In radio-quiet AGN, a'ox 
is likely representative of accretion disc processes, and understand- 
ing the origin and distribution of qJ,^ characteristics will lead to 
a bettering understanding of the primary emission mechanism in 
AGN. 

With the data at hand, we were able to measure the spec- 
tral index of the power law between 2keV and 1090 A (the rest 



frame peak emission observed in the UVW2 filter): aox = 
0.445io(?(/2fcev /fiQggjJ- The spectral slope between 1090 A and 
2500 A is then: q„ = 2.783log(f 'If From these 

^^■'1090 A' ■'2500 A' 

equations we can determine a transformation from our mea- 
sured aox and the standard definition, a'ox, specifically: a^x = 
0.138a„ +0.863Qo:r. 

This leaves us to estimate q„ and thus requires that we make 
some assumptions. Firstly, we assume that the UV continuum 
properties of PG 1407-1-265 are not unusual. While it is clear that the 
emission line properties are unique, McDowell et al. (1995) con- 
cluded that the continuum was rather typical of radio-quiet AGN. 
As such, we adopt the spectral slope of the composite SDSS quasar 
spectrum between 1300-5000 A («„ = -0.44, Vanden Berk et al. 
2001). This brings us to the second assumption, that the slope be- 
tween 1090 A and 1300 A is not drastically different. This is likely 
invalid (e.g. Laor et al. 1995; Shang et al. 2005), but since the sys- 
tematic uncertainty will be introduced to both measurements, the 
differences between the two epochs should remain precise. 

During the low-flux observation (revolution 373), we mea- 
sured a'ox ~ —1.23, which is not atypical of unabsorbed, radio- 
quiet AGN (e.g. Strateva et al. 2005; Yuan et al. 1998; Elvis et al. 
1994). Interestingly, during the high-flux, revolution 206 observa- 
tion, a'ox ~ —1.09, very flat, even in terms of BL Lac spectral 
indices (e.g. Perlman et al. 2005). 

We recalculate a'ox during revolution 206, but this time we 
considered only the 2 keV emission from the hard power law com- 
ponent of the spectrum. This time, a'ox ~ —1.39, much more com- 
parable to the revolution 373 measurement. This may be implying 
that the high energy component originates from a less variable, ac- 
cretion dominated process, whereas the steep, soft component may 
be beamed emission. 



5.2 Long-term X-ray variability 

Considering its high luminosity and radio-quiet nature, the X-ray 
spectrum of PG 1407-1-265 displayed impressive spectral changes 
over an eleven month period. In Figure0some of the X-ray spectral 
difference between the two observational epochs are portrayed. A 
ratio spectrum between the high-flux (rev. 206) and low-flux (rev. 
373) data illustrate how the flux difference between the two epochs 
is not uniform across the spectrum, but gradually increases toward 
lower energies. 

The difference spectrum (high — low) also indicates complex 
spectral changes. A power law (F = 2.18 ± 0.03), corrected for 
Galactic absorption, and fitted to the difference spectrum between 
the rest frame 2.5 — 19.4 keV (1.3 — 10 keV in the observed frame) 
band gives a good result. This high-energy band was selected as it 
appears to show the least amount of spectral variability according 
to the analysis in Sect.|3|and|4| However, there are some deviations 
at the highest energies. On extrapolating the power law to lower 
energies significant deviations are seen (Figure 0, indicating that 
there are changes in the spectral shape. This is supported by the 
different phenomenological fits found when fitting the two data sets 
separately (see Table|2|and|4j. 

The residuals in the difference spectrum imply changes in the 
shape of (at least) one spectral component. The fact that the residu- 
als form an excess at higher and lower energies may suggest that if 
the variations are due to one spectral component, this component is 
likely not smooth in nature (like a power law), but may have some 
curvature (as may be expected from jet emission). 
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Figure 7. In the top panel the 0.25 — lOkeV pn count spectra of 
PG 1407+265 during January 2001 (rev. 206; top curve) and December 
2001 (rev. 373; bottom curve) are displayed. In the middle panel is the ratio 
spectrum of the high-flux (rev. 206) data over the low-flux (rev. 373) data. 
The dashed line marks the average flux change between 0.25 — 10 keV. The 
flux increases more dramatically as energy decreases. The bottom panel dis- 
plays the ratio of the difference spectrum (high - low) fitted with a power 
law(r = 2.18±0.03) in therestframerange 2.5-19.4keV (1.3-lOkeV 
in the observed frame). The difference spectrum has been corrected for 
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Figures. The flux-flux plot for the 0.2 - 0.7keV and 0.7 - lOkeV light 
curves. The higher flux (revolution 206) data can be well approximated by 
a linear (solid line) or a quadratic (dashed curve) fit (shown in detail in the 
inset). On extrapolating the quadratic and linear fits to lower count rates, 
we find that the quadratic model predicts the lower-flux (revolution 373) 
behaviour better. 

5.3 Flux-flux plot analysis 

We further investigated the spectral variability by considering the 
correlation in the flux variations in the two energy bands used. Tay- 
lor et al. (2003) demonstrated that the flux-flux plot can be effec- 
tively used for this purpose. If the varying component of a spec- 
trum has a constant shape (i.e. variations are due to changes in nor- 
malisation only) then the flux-flux relation will be linear. On the 
other hand, the relation will not be linear if the spectral shape also 
changes. 

We constructed a flux-flux plot (Figure|8} as demonstrated by 
Taylor et al. (2003). In considering the variable high-flux observa- 
tion, both a linear and quadratic function fit the data equally well 
(X^ ~ 0.93 and 0.37, respectively), although the quadratic fit does 
seem a better approximation to the data (see the inset in Figure[8). 

The possible non-linearity displayed in the flux-flux plot dur- 
ing revolution 206 implies that the short-term (within the observa- 
tion) variations are not due to a simple change in normalisation of 
one component, but require a change in the spectral shape as well 
(i.e. pivoting of the power law at high energies; e.g. Zdziarski et al. 
2002). This is consistent with the flux-dependent spectral variabil- 
ity seen in Figure|4|and the spectral differences seen in the analyses 
of the high- and low-flux state spectra during revolution 206 (see 
Figure|5| Table|3| and Section l3.2.3> . 

Extrapolation of the linear and quadratic fits to lower count 
rates in the flux-flux plot showed that the quadratic fit also predicted 
the low-flux state of revolution 373 much better than the linear fit 
(Figure|8}. As with the difference spectrum (Figure^, this further 
emphasises that the long-term variability, which is only examined 
at two epochs separated by eleven months, required changes in the 
shape of the variable component. 



6 DISCUSSION 
6.1 General Findings 

The main results of this analysis are the following: 
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(1) In comparing X-ray data from two epochs, PG 1407+265 ap- 
peared to exhibit two distinct states: a variable, high-flux state (rev- 
olution 206); and a quiescent low-flux one (revolution 373). At both 
epochs, a power law component with F ~ 2.24 was required, while 
during the high-flux state an additional, steeper (F ~ 2.5) power 
law was also needed. Significant spectral and flux variability was 
seen during the high-flux state, whereas during revolution 373, the 
variability was negligible. 

(2) Short-term variability (within revolution 206) and long-term 
variability (between the two epochs) was consistent with spectral 
hardening with diminishing intensity. In comparing the hard (0.7 — 
10 keV) and soft (0.2 — 0.7 keV) light curves during the high-flux 
state, there appeared to be a short lag such that the hard energy band 
followed the soft by about 3000 s. 

(3) The radiative efficiency measurement during the bright, rev- 
olution 206 observations was r) > 0.29. The value is consistent 
with emission from an environment around a Kerr black hole or 
relativistic beaming. 

(4) Though not statistically significant, there appeared to be some 
UV variability during the high state, which showed similar trends 
as the X-ray light curve. 

(5) The optical-to-X-ray spectral slope was estimated at both 
epochs from simultaneous UV data, and found to be significantly 
flatter during the brighter X-ray phase. UV variability between the 
two epochs was within about 10% and typical for quasars, whereas 
the average broadband X-ray flux changed by about 60%. 

6.2 The physical nature of the X-ray emission 

The broken power law and blackbody plus power law models pre- 
sented thus far are informative. For example, both models illus- 
trated significant spectral variability between the two epochs, as 
well as on shorter time scales during revolution 206, but they lack 
physical intuition. A sudden break in a power law, or a high and 
variable blackbody temperature, appears unphysical. 

In this section we attempt to derive a physically motivated 
model for the spectral behaviour of PG 1407+265. In doing so, 
we demand that all the different spectra can be fitted by a physi- 
cally consistent model. This means that we not only attempt to de- 
scribe the spectra at both epochs with a consistent model, but also 
the high- and low-flux spectra during revolution 206. Consequently, 
the fits in this section include three rather than two data sets (high- 
flux revolution 206; low-flux revolution 206; lowest flux revolution 
373). For simplicity, we use only the pn data in this section. 

As opposed to a broken power law, a double power law does 
have physical foundation. Initially we fit all three spectra with a 
double power law model. We required that the photon indices were 
constant in each spectrum and allowed only the normalisations to 
vary. This provided a good fit to the data (x^/dof = 1.00/1607) 
and illustrated how the hard power law component dominated dur- 
ing the lowest flux state (revolution 373) and the steeper component 
became more relevant with increasing flux. 

We continued to build on this by also fixing the normalisation 
of the hard power law to be constant in all three spectra. There- 
fore, the only variable component at each epoch was the normal- 
isation of the soft power law. This provided a slightly worse, but 
stifl acceptable fit (x^/dof = 1.04/1609). The hard power law 
component (F — 2.24:'^'^^'^) was essentially established by the 
revolution 373 spectrum. The contribution of the steep power law 
(F = 2.55 ± 0.02) component to the total spectrum during revo- 
lution 373 was negligible. However, as the flux increased, the con- 
tribution of soft component became more important. In the highest 



flux state (high state during revolution 206) the steep component 
was attributed to 62 ± 1% of the unabsorbed, 0.25 - 10 keV flux, 
diminishing slightly to 58 ± 1% in the low-flux spectrum of revo- 
lution 206, and completely negligible (0%) at revolution 373. 

It is interesting that this simple double power law model pro- 
vides a reasonable fit to three different spectra. As seen through- 
out the paper, the soft component was more variable on short and 
long time scales. The fastest time scales at which we would ex- 
pect to see variations associated with the accretion disc, would 
be on dynamical time scales {td = 2 x IQ^ Ms{^)'^/'^ s). For 

PG 1407+265, with a black hole mass of ~ 10^ M© (Mg = 10) 
and r = Tg = GM/(? this corresponds to variation on the order 
of 20 ks in the observed frame. We are clearly seeing variations on 
much shorter time scales. 

Emission from an accretion disc could be consistent with 
the harder power law component. The absence of clear signatures 
from the disc (e.g. iron emission line, reflection component above 
lOkeV, and variability) is not so surprising given the high lu- 
minosity (> 10**^ erg s"^ ) of PG 1407+265 (e.g. see Reeves & 
Turner 2000). In fact, the harder power law in the double power 
law model can be successfully replaced with a thermal Compton- 
isation (Titarchuk & Mastichiadis 1994) component (x^/Aoi = 
1.05/1607). In this case, the steep power law remained unchanged 
compared to the double power law fit. The seed photons have a tem- 
perature of kT ~ 33 eV and the plasma has a temperature of about 
56keV. The situation is not inconsistent with predicted physical 
conditions, although the measured model parameters are not well 
constrained by the fit. 

It is conceivable that the soft, variable power law component 
can be attributed to emission from an inner jet. The slope of the 
X-ray spectrum (F ~ 2.5) is steeper than the radio spectrum (e.g 
Ulvestad et al. 2005), which argues against inverse Compton scat- 
tering as the emitting process. The steep spectrum and remarkable 
variability strongly favour synchrotron radiation (e.g. Wilson & 
Yang 2002). 

The jet could be working intermittently: present and dominant 
during the high-flux revolution 206 observation, and absent during 
the later, low-flux observation. In the spectral fits presented in this 
section, we have demanded that the spectral slope remains constant 
at all epochs. This is, of course, not a necessary constraint as the 
X-ray slopes of jet emitting AGN are know to vary dramatically 
on various time scales. In fact, based on spectral modelling, we do 
not rule out variations in the photon index of the jet component 
between epochs, but note that the lack of variability seen during 
the low-flux state, makes the argument pointless. It should also be 
noted that the jet emission process is probably better described by 
a curved spectrum in the XMM-Newton energy range (Perlman et 
al. 2005). Variations in this curvature can be responsible for some 
of the spectral variation implied by the difference spectrum and the 
flux-flux plot analysis (Figure0and|8| respectively). 

Combined jet and accretion disc models have been fitted to a 
number of nearby, radio-loud, unbeamed objects (e.g. Cen A, Evans 
et al. 2004; NGC 4261, Zezas et al. 2005), and have also been pro- 
posed to account for the apparently higher X-ray flux in beamed 
objects (e.g. Browne & Murphy 1987; Kembhavi 1993). The fact 
that neither component in PG 1407+265 appears to suffer from ob- 
scuration, favours the conception that PG 1407+265 is beamed as 
well (consistent with its radio jet). If the presumable lag detected 
between the two components (Figure|3} is attributed to light travel 
time delays, then the two components (likely jet base and accre- 
tion disc corona) are separated by only about 1 r-g. This essentially 
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places the base of the jet inside the corona, giving credence to the 
idea that a jet could be replenishing the plasma reservoir with high- 
energy particles. 

6.3 The spectral energy distribution 

The spectral and short-term timing behaviour do support the notion 
that PG 1407-1-265 is being viewed pole-on (Blundell et al. 2003). 
This can, in principle, describe the weakness of the high-ionisation 
emission lines seen in the optical and UV if, as Blundell et al. de- 
scribe, we are looking directly into the broad line region. On the 
other hand, the strong, low-ionisation Fe II emission could be colli- 
sionally produced at large distances where the radio jet components 
are found (~ 100 pc from the core, Blundell et al. ). 

The suggestion that the X-ray jet works intermittently also de- 
scribes why the SED is consistent with an unbeamed continuum, 
assuming that the jet is more often off than on. This follows from 
the dramatic differences in aox, as measured from simultaneous X- 
ray and UV data, which indicate that the continuum can certainly 
be much flatter (apparently beamed) when the jet is on. In addition, 
if the X-rays are assumed to originate from the base of the jet and 
hence feeding the radio jet, the irregularity of the X-ray jet emis- 
sion could help explain the weakness of the radio emission relative 
to the optical. 



7 CONCLUSIONS 

We have presented two XMM-Newton observations, separated by 
eleven months, of the radio-quiet quasar PG 1407-1-265. The X-ray 
behaviour is consistent as arising from a combination of a rela- 
tivistic jet and accretion disc processes. The spectral and timing 
variability strongly suggest that the jet works intermittently, which 
could help describe some of the complexities seen in the SED of 
the quasar. The variations found in aox, measured from contem- 
poraneous X-ray and UV data, demand the need for simultane- 
ous, multiwavelength observations to accurately define the SED of 
PG 1407-^265. 
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